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Abstract. In the framework of the Bogolyubov approximation and using the 
Bogolyubov inequahties we give a simple proof of the coexistence of two non- 
conventional Bose-Einstein condensates in the case of some superstable Bose sys- 
tem whose atoms have an internal two-level structure and their energy operators 
in the second quantized form depend on the number operators only. 



1. Introduction 

The phenomenon of Bose-Einstein condensation (BEC) was described first by 
Einstein pQ. Until recently the best experimental evidence that BEC could occur in 
a real physical system was the phenomenon of superfluidity in liquid helium as sug- 
gested originally by London who introduced the concept of macroscopic occupation 
of the ground state [2] and conjectured that the momentum-space condensation of 
bosons is enhanced by spatial repulsion between the particles [3]. However, nowa- 
days, there exists a considerable amount of experimental evidence for BEC [U |5]. 

Given the difficulty of the problem of proving the existence of BEC from a math- 
ematical point of view, it is desirable to have idealized models in which one can 
develop concrete scenarios for BEC. In this sense, we are interested in studying 
Bose systems whose energy operators consider repulsive mean interactions repre- 
sented by diagonal operators in the occupation numbers. It frequently leads to 
thermodynamically stable systems which can be classically understood. 

We are considering a class of systems for which various Bose condensates, in the 
sense of macroscopic occupation of the ground state, coexist |6l[7]. Specifically we 
are interested in the study of thermodynamic behavior of some systems for which 
the Bose atoms have an internal two-level structure. We prove that the existence of 
two ground state levels, one of them with negative energy, leads to an enhancement 
of condensation. 

This paper is divided as follows. In section [2] we describe the basic mathematical 
notions associated to this kind of systems and we study the stability of a system of 
atoms with internal structure. In section [3] making use of the so-called Bogolyubov 
approximation [HI [9] we give an extremely simple proof of non-conventional BEC in 
the sense of macroscopic occupation of the ground state [TOl [TTl [12] . We prove that 
the existence of two ground state levels, one of them with negative energy, leads 
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to an enhancement of condensation in the sense of macroscopic occupation of both 
internal states [HI E] . The proposed approximation (Bogolyubov approximation) is 
obtained by substitution of the zero mode creation and annihilation Bose operators 
by suitable chosen c-numbers. In this proof a significative role play a well known 
variational Bogolyubov inequalities and the fact that the involved Hamiltonians 
are written in terms of the number operators. It enables us to give a simplified 
demonstration of thermodynamic equivalence of the limit grand canonical pressures 
corresponding to the energy operator of the system and the respective approximating 
Hamiltonian. 

Atoms BEC's in two different hyperfine states confined in a single trap with a 
time- varying Raman coupling between the two levels are studied in Ref.l3. The 
Hamiltonian associated to such system, being a, b the annihilation Bose operators 
corresponding to the two different states and under the two mode approximation is 
given as: 

H = Ha + Hh + Hint + -f^las, (1) 
Ha = c^aci^d + -^d^a^aa, (2) 

Hb = toM + ^PPhh, (3) 

Ant = ^atafetfo, (4) 

H\^ = fi(t)(e***a^6 + e~'H^a), (5) 

where Ha, Hb describe systems with self interactions. Hab is associated to collisions 
between both systems and Ua, Ub, Uab are constants. Finally H\as represents the 
operator associated to the Raman coupling. 

We shall prove by using techniques of Bogolyubov approximation that an anal- 
ogous scenario for the coexistence of different condensates can be obtained in the 
case of superstable systems of two-level Bose-atoms, described by energy operators 
of the mean field type, undergoing self interactions and collisions. Obviously the 
Raman coupling is excluded from the respective Hamiltonians. 

2. The Model 

The one-particle free Hamiltonian corresponds to the operator S*' = — -y defined 

on a dense subset of the Hilbert space H} = L^(A/), being A/ = [— |, ^j"^ C M"' a 
cubic box of boundary dAi and volume Vi = l'^. In other words, the particles are 
confined to bounded regions. We assume periodic boundary conditions under which 
S*' becomes a self-adjoint operator. 

In this section it is assumed that the Bose atoms have an internal two-level struc- 
ture analogous to the SU(2) spin symmetry [U [7]. In this case any one-particle 
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wave function has the form ® s where, G L'^i^i) and s G represents the 
internal state. Therefore the vector space associated to this system is in fact, 

We shall study the model of Bose particles whose Hamiltonian is given as: 



Hi = H? + U! 



H? + ^° 



it 



7 




(6) 



where g[x) > ax'^,a > is a continuous function and cr = + or — depending 
on the corresponding label of internal energy. The second term at the right hand 
side of eq. ([6]) represents the intrastate collisions, or the self-scattering, the third 
term represents the interstate collisions or the cross scattering. The sum in ([6]) 
runs over the set A^* = {p G M'^ : = 27™^//,^^ = 0, ±1, ±2, , ...a = l,2,..d}, 
cip cr are the Bose operators of creation and annihilation of particles defined 



p,cri 



on the Fock space J-'b and satisfying the usual commutation rules: [d, 



^q,0'l^p,(T2 ^p,(T2^'l'<^l 



q,o-u "'p,(T2 



'^p,q'^(Ti ,(72 • ^PjC 



d' a-^p,o- is the number operator associated 



to index p and internal label a and 



p = 0, 
Pt^O. 



(7) 



In this case Hi 



total number operator and 
with exclusion of no cr. 



Ea.peA* ^«,a(p)^p,<7 and 70,7 > 0. iV = ^ 



at „dr. rr IS thc 



,peA*\{o} "p.o-^p 



cr.pGA,* "-p.o-^-p.o- 
,dp^o- is the total number operator 



Note that the boson Fock space J-'b is isomorphic to the tensor product (8>cr,pGAf 



p,cr 



where J^^^^ is the boson Fock space constructed on the one-dimensional Hilbert space 



n 



P,<T - 

Let 



{Oe 



ip X 



ecrjeGCj where e_ = (0, 1) and e+ = (1,0). 



PiiP, fx) = — In Tr^3 exp 



-l3{Hi - fxN) 



be the grand-canonical pressure corresponding to Hi, where P = 9 is the inverse 
temperature. If Hi{ij) = Hi — fxN, the equilibrium Gibbs state (grand canonical 
ensemble) (—)^ , . is defined as 



(A) 



Tr^s exp [-l3Hi{fi) 



-1 



Tr 



Aexp{-l3Hi{fi) 



(9) 



for any operator A acting on the symmetric Bose-Fock space. The total density of 
particles p(/i) for infinite volume is defined as 




Vi 



lim pi{fi) = p{fj), 



(10) 
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and the density of particles Pp^aiP, Ij) associated to the energy label Acr(p) is defined 
as 



J™ ("S^/ = Jim Pp,<x,/(/^>/^) = Pp,<x(/5,/^)- (11) 



We shall say that the system undergoes a macroscopic occupation of the single 
(p, o")-mode particle level if pp,o-(/5, yu) > 0. 

In the case of models with a two label internal structure, [6], [7] two macroscopic 
occupations independent on temperature (non- conventional BEC) in the sense that 



Po,-(/^)>0, Po,+ (p)>0 (12) 



can coexist. 

We are interested in proving the superstability of the model, in other words we 
verify that the limit pressure lim pi{f3, p) exists for all p G M. A useful criterion [Ti] 

says that we are in presence of a superstable system if given an interaction f// defined 
on the Bose Fock space J-'b, there exist two constants Ci > and C2 > such that 



Ul>-^N + ^N\ (13) 



Then, the superstability of the system whose energy operator has the form is 
ensured by the following result: 

Proposition 2.1. The operator Uj given by Ij^ satisfies the inequality 



Proof. For x, y G M we have x"^ + y"^ > (x + yY /2. Moreover, since (0'p,CT)^(ap,o-)^ = 



1) > and the spectrum Sp(np^o-) = N U {0}, for all a = ±, p G A^*, we 



obtain, in the sense of operators, the following sequence of inequalities: 
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7o^> , 70.2 , 70.2 , A^'f 



>-j-N + j-nl_ + j-hl^ + a 



Vi Vi Vi "'^ Vi 

- Vi 2Vi 2V^ L 0,+ ^ J 

^_2o^^^mn|a^^,_ (15) 
This concludes the proof. □ 



3. Approximating Hamiltonian 
For any complex number c e C, let consider the coherent vector 0c,ct in ^o,a, given 

as 

oo ^ 

0c,. = e-^'HV^^-L(^ic) (ag"0o,., (16) 

n=0 

where 0o,(t = 0o ® is the vacuum vector in J-'o,a- Then, for any operator A defined 
in J^B, it is possible to construct an operator A{c) such that 

h[,A{c)%) =(^;®0e,a,M®0c,a) , (17) 

where ^2 G being J^'^ — '^a,pGA*\{o}^p,a- The transition from operator A to 
operator A (c) is known as Bogolyubov approximation and it consists in replacing the 
operators Oq and Oq in any A expressed in the normal form ( all creation operators 
are grouped on the left hand side of the group of annihilation operators in A) 
operators with the complex numbers and \^ c. Then, in the framework of the 
Bogolyubov approximation, replacing the operators aj,^ and ao,cr in Hi{fJ^) with the 
complex numbers \/Vic, \/Vic, for a — — and \/Vif}, \/Vir), for a = +, respectively, 
we get the following approximating Hamiltonian H^^^^{c,r], /i), 
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a,pGA,*\{0} '■ a,peA*\{0} 



-Vig (f ) + (\_(0) - fi)Vi \c\' + 70^ \c\ 



+ (A,,+ (0) - fi)Vi \r]\' + ^oVi + iVi |c|' \r^\^ , (18) 

with c, G C. 

Note that Hf''^^'^ {c, fi) can be rewritten in a more suggestive form as: 

' a,peA*\{o} 

+ (A;,+ (0) - fi)Vi \vf + 70^ 1^1' + jVi \cf \r]f , (19) 



where 



H'^'^'ifi)^ Y Mp)-f^X^ + Vig{^] (20) 



a,pe\*\{0} 



denotes the so-called mean field Hamiltonian with excluded zero mode. 
In this case we have the following theorem 

Theorem 3.1. The following inequality takes place: 

Tr^. e-^^r^'^^'"''^) < Tr^, e'^^'^^). (21) 
Proof. From the above definitions it follows that 

Tr^^ (^e-f'Mcr^,!^)^ = sup Y e^P (^;, Hi (c, r/, /i) 

= sup Vexp f-/5 ® ® 0^,-,^z(/i)^; ® 0c,- ® 



< Tr^Je-^^'(^) ) . (22) 



The supremum is over all orthonormal basis of the corresponding Fock space JF^. 
The first inequality is the Peierls inequality; the second one is obtained taking 
the coherent vector 0c,- ® 0?^,+ as one of the vectors in the orthonormal basis of 
^0 = -^0,- ® -^0,+- This complets the proof. □ 




(23) 
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3.1. Bogolyubov Inequalities. Let Ha,i and Hb^i be self-adjoined operators de- 
fined on P C J-'b- Pa,i{P,fJ'), Pb,i{P,fJ') represent the grand canonical pressures 
corresponding to the operators Ha^i, Hb^i- In this case the following well known 
Bogolyubov inequalities 

/ ,,\\ 

V, 

hold, where (— ) ^^^-j, (— ) jj-j^ j(^) are the Gibbs states in the grand canonical ensemble 
associated to the Hamiltonians Ha^i,Hf,j, respectively. 

3.2. Global and Local Minimums for Convex Functions. 

Proposition 3.2. Let G : U ^ W be convex on a convex set U G L ( L is a normed 
linear space). If G has a local minimum at x, then G{x) is also a global minimizer. 
The set V ( conceivably empty ) on which G attains its minimum is convex. And if 
G is strictly convex in a neighbourhood of a minimum point x, then V = {x}; that 
is , the minimum is unique. 

The proof is contained in Ref.l5. 

3.3. Equivalence of Limit Grand Canonical Pressures. 

Theorem 3.3. Under any of the conditions 

(i) 270 > 7, > A(27o + 7)/(27o - 7), 

(ii) 7 = 0, G M, 
(Hi) 270 = 7, /i e M, 

the systems with Hamiltonians HiifJ^) and H^^"^^ {c, 77, /i) are thermodynamically equiv- 
alent in the following sense, 

lim sup pfPP''(/5,c,r/,/i) =p(/?,/i), (24) 

|c|,|r)|:c,»;GC 

where pf^\l3, c, 77, /x) = {(3Vi)-^ InTr^^ exp -/5(^f''P'(c, r/, /i)) 

Proof. Applying the Bogolyubov inequalities for pressures (123!) and inequality (1211) . 
we obtain that: 

< pKA/i) -Pr'(/?,c,r/,/i) < i(^fPP^(c,r/,/i) - = ^Kc,r/,/i). (25) 

The right hand side of inequality (l25l) can be rewritten as, 
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Mc T], /i) = (\_(0) - fi) \cf + 70 \cf + (A;,+ (0) - fi) \r]f + 70 l^^l' + 7 \cf \vf 



lPo,-Po,+ ) (26) 



where po,a = ^o.^/VJ. 

(i) Under the condition 270 > 7 and according proposition 13.21 g(x. y) : M."^ —>■ M. 
given as: 

gix, y) = (A/ _(0) - + 7ox^ + (A/,+(0) - /i)?/ + 70?/^ + 7X?/, (27) 
is a strictly convex function attaining a global minimum at 

47o - 7' 



^0 7-2 — > l^»J 



_ 7(A^,-(0)-^)-27o(A,,+(0)-/^) 

^:;t^2 ' (29) 

with {x,y) G M^, since d^g = 270 > and d^gdyQ — [dlyg^ = ^ll — 7^- Taking 
Ai__(0) = —A and Ai_+(0) = A, we have that g{x, y) = —(A + /x)x + 7oX^ + (A — fi)y + 
lov"^ + l^y and 

^o = \co\' = :^ + :^, (30) 
270 - 7 270 + 7 

yo = \Vo\' = -:^^ + :^^, (31) 
270 - 7 270 + 7 

for /i > A(27o + 7)/(27o — 7). It is easy to verify that: 

gixo, yo) = - (tt^ + TT^. — 1 • (32) 

V270 - 7 270 + 7/ 

On the other hand, being /i; = /i + 7o/VJ we have 

(^r(Po,-,Po,+))H,(^) = ( ( 5^(A/,<x(0) -^OPo,a + 7oPo,a) + lPo,-Po,+ ) (33) 
where 

(7;(x, y) = (A,,_(0) - ^ll)x + 7oa;2 ^ (Az,+(0) - /^Oz/ + 7oy' + 72:2/. (34) 

With Xi-{0) = — A, Ai,+(0) = A, and under the same before estabhshed conditions 
gl has a global minimum at the positive real values 

^0, = + (35) 

270 - 7 270 + 7 

* A /i; 

yo,i = -7> + 7^ — — • 36 

270 - 7 270 + 7 
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Clearly Xqi ^ xq and yli yo when VJ -> oo and lim gi{xli,yli) = g{xo,yo)- 
Moreover, gi{x,y) < g{x,y) for all {x,y) E M+^, then 

inf Q*i{x,y)< inf g{x,y) (37) 



which implies that 



inf g{x,y)+ sup {-^;(x, y)} > 0. (3J 



With these definitions we can rewrite Ai{c,ri,fi) as, 

Az(c,r/,/i) = ^(|r/|% |c|^) + {-g!{po,-, Po,+)) H,i^^)■ (39) 

The energy operator ([6]) is diagonal respect to the number operators. Therefore, 
since the spectra of these operators coincide with the set of non negative integers, 
this model can be classically understood by using non negative random variables 
defined on a suitable probability space Qi. Let Qi be the countable set of sequences 
uj = {uj{p, a) eN : a = ±,p E AJ} C N U {0} satisfying 

J2 ^(P,^) < oo . (40) 
o-,peA* 

The basic random variables are the occupation numbers {?T-p,cr '■ j = 1, 2, a = ±}. 
They are defined as the functions Up^a : fli ^ N given as np^„{u) = uj{p, a) for any 
UJ E Qi. The total number of particles in the configuration u is denoted as N{uj). 
Then the total number, excluded the zero mode is denoted as N'{uj). 
Let Hi{fi) be the function of the random variables Up^a defined as 

Hiip) = I E (^'.'^(P) - /^')^P.'^ + y^<'^ ] + ^^o,-r.o,+ + Vwi^). (41) 

\a,p6A* I J I I 

Let P be a probability density defined for any uj E Qi as 



UJ 



n -1 



5]exp(-/?[/7,(/i)]H) 



exp{~P[Hiifx)]iuj)). (42) 



Being E the expectation related to P, the concavity of —gi implies that this 
function has a global maximum, then we obtain. 



(-^r(Po,-,Po,+))//;(^) = E(-^;(po,-,Po,+)) < E( sup {-gi{x,y)}), (43) 

(x,y)GK+2 
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where pco- = "^0,0-/^/. From this it follows that 

< p/ /u) - sup p^PP' {(3,c,rj,ii) < inf A/ (c, 7;, /x) 

|c|,|?7|:c,r?eC |c|,|?j|:c,?7eC 

-,M^/^f ^^'(|c|^ 1^/1^) + sup {-g;{x,y)}. 

\c\,\ri\:c,rieC {x,y)eR+^ 

Therefore, 

< lim inf Ai(c,r/,/x) < g{xo,yo) + lim {-^r(4,^ l/o,i)} = 0- 

Vj— >oo |c|,|r;|:c,??eC Vj— >oo 

(ii) We have that: 

V, ^) - Hi{^))h^^^~^ = -(A + /x) |c|^ + 70 |c|^ 
+(A - IJ) + 70 + ((A + /i)po - - 70P0 -(po - - 

+((// - A)po,+ - 7oPo,+(Po,+ - 

Let po,(T = {Po,a)Hi(p)- From the above result and noting that 

inf Api{(3, IJ) - pf^\(5, c, r],iJ,)}^ pi{(3, /x) - sup pt^^\P, c, //) 

|c|,|r?|:c,r/eC |c|,|?j|:c,r;GC 

and {{po,ay)H,(f,) > {{po,a) H^^^)f it follows that, 

< pi{(3, //) - sup p^^^'{(3, c, r],ij)< inf {-(A + //) |c|^ + 70 |c|^}+ 

|c|;ceC |c|:ceC 

+ inf {(/. - A) Ir;!^ + 70 + (// + A + ^)po,- - 7oP^,-+ 

+(/^ - A + ^)Po,+ - 7oPo,+- 

Note that 

0, /^<(jA-y^, 
\ ^ 

^ - (jA + ) 

/i>(7A-— . 

470 



sup {(/i - (tA + :^)po,a " 7oPo,a} = ' 



PO,cr>0 V'i 



On the other hand 



inf {(crA - + 7oX^} = < ia-aX) 
x>o ^ — '—^ ii> aX 



0, /J, < crA, 

— (7/ 
470 
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The respective maximum and minimum are attained in the following way, 
0, 

fj, — (T\ 
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P0,a 



H<aX- — 
Vi 



, ix>aX-—, 
270 Vi 



0, II < aX 

— aX 



270 



-, fj, > aX. 



(51) 



Combining these results we obtain, in the thermodynamic limit, the proof for the 

case 7 = 0. 

(iii) In the case 270 = 7 we have g{x, y) = —{X + ix)x + 70(0; + yY + (A — ii)y. 
(a) ji e (—00, —A]. In this case g{x,y) > 0, for all {x,y) e (K"*")^, 
then inf {g{x, y)} — g{^-i 0) = 0, where the minimum is attained at 



(x,2/)G(IR+)2 



(b) 11 e (-A, A]. We have. 



xo = 0, I/O = 0. 



(52) 



, 'V'L^M^^y)} ^ , inf ^ |-(A + /i)x + 7o(x + y)2 + (A-//)y} 



inf {gix, 0)} = inf {-(A + /x)x + ^qx"^} 

a;eIR+ a;eIR+ 



470 



where the minimizer is given by, 

/^ + A „ 
xo = — — , yo = 0. 
270 

(c) 11 e (A, 00). Let ri = A + r2 = A* — A. Then, 

9i.x, y) = -rix + 70(2; + yf - r2y. 

In this case, 

inf \q(x, y)} — —. 



The minimizer is given by. 



2^0 = — , 2/0 = 0. 
7o 



This ends the proof. 

Corollary 3.4. Assume that supj.>o {(// — Oi)r — g{r)} exists. Then, 



(53) 

(54) 

(55) 

(56) 

(57) 
□ 



p(/3,/^) 



A^ 



+ 



270 - 7 270 + 7 
/ii(A,/i)+p^^'(A/x), 
L /i2(A,/i)+p^^'(/3,/x), 



+ (/?,//), 27o>7>0,/i> 



A (270 + 7) 
270-7 



7 = 0, 70 7^ 0, 
270 = 7, 



(58) 
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r 0, 



(/i + A)^ 

470 
//^ + A^ 

270 ^ 



H e (—00, —A] 
/X e (A, 00) 



r 0, 



470 
A/x 

To ' 



/X e (—00, —A] 
A* e (-A,A] 

jji e (A, 00). 



p^^' /.) = sup inf { (/. - a)r - ^(r) + p"^ a) | , 

r->0 i^SU 1^ J 



(59) 
(60) 



being p^"^ (/3, a) the limit grand canonical pressure for the ideal Bose gas with excluded 
zero mode at chemical potential o; < 0. 

Proof. Let i^V) = H^^\^A + 'i E.,peAr\{o}(^kJ'<. = H^^' ip) + fD and H^' = 
SCTpeA*\{o} Ai,(7(p)^p^CT^p,o-- Let f[{Q) be the canonical energy at finite volume Vi and 
density g associated with the Hamiltonian F, defined as, 



f[{Q) = -TjlnlV (B,jv) e 



/3r|^(B,jv) 



(61) 



where Ti^^'^^ is the Hilbert space representing a //-particles Bose system. Let Qq 

the density of particles in the zero mode. Let fiig — Qo), ff^^ {q ~ Qo) , fi'^ {q — Qo) be 

the finite canonical energies associated to the operators H',H^^ ,Hi' respectively 

and f {g — go), f^^ {g — go), f^'^ {g — go) their respective limit canonical energies. 

We only sketch the proof which can be found in Refs.11,12 in the case g{x) — ax"^. 
In this case the following identity can be easily verified: 



/ 1 —^lOf) 

fi(Q-Qo)^g{Q-Qo)-^He ) hUs-,o) 



(62) 



where {')Hoi{g-eo)^ represents the canonical Gibbs state at density g — go associated 
to the ideal gas with the single mode zero excluded. Using the following facts, 



(e , 



'I <1, 



and noting that 




= 0, 



(63) 
(64) 

(65) 



H? {8-00) 
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we get 

lim -577 In e =0. (66) 

Vi->oo pVi ' H'^ [Q-eo) 

This implies that in the thermodynamic hmit 



Then, 



but 



Therefore, 



f{Q- Qo) = n {0 - Qo) = g{Q - Qo) + r {q " Qo)- (67) 



P iP, /i) = sup {fi{g - ^0) - /^^ } 
snp{fi{g - Qo) - g{Q - Qo) - f"^ {q - Qo)}, (6J 



f ^ {q - Qo) = snp{a{Q - qo) - p'" (a)}. (69) 

a<0 



p'{(3,fi) = sup sup{(/i -a){Q- Qo) - g{Q - Qo) + p'"^ («)} 

= supsup{(;U — a){r) — g{r) + p^'^ («)}• (70) 

r>0 a<0 

On the other hand it is clear that in the case 270 > 7 we obtain, 

p(/?,/i) = sup {(/i + A)x - 7oX^ + (/i - A)y -7o?/ - 72;y} +p'(/5,/i) 

(x,y)G(M+)2 

= (yu + A)xo - 702^0 + (/^ - -^)yo - loVo - l^oVo + P (/3, /x) 

^^ + 7/-—)+P'((^^^')- (71) 
270-7 270 + 7/ 

In the case 270 = 7 the proof follows from the theorem I3.3[ 

□ 

3.4. BEC. 

3.4.1. Coexistence of two non- conventional Bose-Einstein condensates. In this sec- 
tion we theoretically predict the exciting possibility of coexisting non-conventional 
Bose-Einstein condensates. 

Theorem 3.5. Under the conditions 270 > 7 and fi > X{2jo + 7)/(27o — 7), the 
system given by ^ displays simultaneous non- conventional Bose-Einstein conden- 
sation of the two levels associated to the zero mode and the amounts of condensate 
are given as: 

PoAp) = + (72) 

270 - 7 270 + 7 
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Po,+ (/i) - 



A u 

+ 



270 - 7 270 + 7 



Proof. p{P, n) can be rewritten as, 

1 / ia- — a+Y 



4 V 270-7 



(73) 



(74) 



270 + 7 

where a_ = A + /i, = /i — A. Using the convexity of Pi(/3, /i) and pf^'''^(/3, c, r], /x) 
respect to a_ and a+ we get from the Griffiths theorem [16] and theorem 13.31 



Po,<7 = ,lini da^pi{(3, fi) = da^ hm pi{(3, fi) = d^^ hm sup pi{(3, c, 77, n) 



1 / (a- — («_ + 0;+)^ 



+ 



(75) 



^4 V 270-7 270 + 7 

for cr = — , +. An analogous procedure leads to eq. (173|) This yields to the proof of 
the theorem. □ 



Corollary 3.6. For 7 = we have, 

0, yU G ( — C)0, —A] 



/i + A 

270 ' 



Po, 



and for 7 = 270, 



Po, 



/i G (-A,oo), 

/i G (—00, —A] 



jl — \ 



/i G (— C)0, A] 
, /iG(A,cx)) 



(76) 



r 0, 

p + A 

270 
P_ 

7o' 



/i G (-A, A], 
P e (A, 00], 



Po,+ = 0, p G 



(77) 



Proof. This proof is a direct consequence of Corollary 13.41 



□ 



Let Tss^i{l3, fi) and Tcs^i be thermal averages related to the self-scattering and cross- 
scattering operators defined as: 

^ss,/(/5,p) = V;~^ ^{{ai^aYalJfj^^^-^, T^s,iW,fi) = (po,+Po,-)h;(^)- (78) 
Corollary 3.7. For 270 > 7 and jj, > A(27o + 7)/(27o - 7), 



lim Vi Vcs,/(/3,p) = lim {po,+Po,^) = Jim (Po,+)j^,(^)(Po,-)i^,(^) 
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lim y-\,//3,/i) = 2( (-A 



and for 7 = 0, 



lim V,s,/(/?,/x) = < 



7 



270 
I 270^ ' 



/i 



270 + 7 



/i G (—00, —A] 
fi E (-A,A] 

/i G (A, 00). 



Corollary 3.9. For 7 = 270 the following identities 



lim V;-^(rss,z(/?,/i) + 2r,,,i(/?,/i)) = < 



+ A 



A/i 

2"' 
I 7o 



/X G (—00, —A] 
^ G (-A,A] 

H G (A, 00) 



hold. 



(81) 



f82) 



4. Generalized BEC 

We consider the case g{x) = ax'^, a > 0. In this case Hamiltonians ([6]) and Hf^^ (/i) 
become 



peA;\{0} 



V, 



Vi 



(83) 
(84) 



a,peA*\{0} ^ 

respectively. 

Let p'liP, /i), pf^^ (P, fj) be the grand canonical finite pressures associated to these 
operators. For every /x G M, in the thermodynamic limit, it is proved in Ref.l7 that 



lim p'iiP,fx) =p'(/?,/i) = lim pj^^ 
= p^^' (/?, /i) = sup L^'iP, a) + ^^^^ 



^5) 
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Proposition 4.1. Let 



(2vr) 



[0,00)" 



I l|2 

IpII 

2 -1 



(86) 



and 270 > 7, then 



, , MF'm N A(27o + 7) ^ a 



270-7 270 + 7 270-7 2 



270-7 a 270 + 7 270-7 

(87) 



Moreover this model displays absence of macroscopic occupation of nonzero single 
levels: 



lim/^; =0 (88) 
for 0" = +, — and p G A;*\{0}, but generalized BEC can be verified in the sense that 



lim lim y ( ^ ) > 0, (89) 

a,peA*:\i{p)<S ^ ' ' 



being fip = fip^ + hp^^ 



5. Conclusions 

The pressure of the model given by ([6]) can be exactly determined for a range of val- 
ues of the chemical potential /i by using the energy operator f|T9l) in the framework of 
the Bogolyubov approximation. The superstable system displays non-conventional 
BEC consisting in the macroscopic occupation of both ground state levels. The con- 
densates can eventually coexist depending on the values of fi. Moreover, as a subtle 
consequence of these results we have determined the exact values of some thermal 
averages associated to the cross scattering and self-scattering operators. Finally, in 
the case g{x) = ax"^, a > 0, generalized BEC holds. 
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